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We analyze the quantum walk in higher spatial dimensions and compare classical and quantum 
spreading as a function of time. Tensor products of Hadamard transformations and the discrete 
Fourier transform arise as natural extensions of the "quantum coin toss" in the one-dimensional walk 
simulation, and other illustrative transformations are also investigated. We find that entanglement 
between the dimensions serves to reduce the rate of spread of the quantum walk. The classical limit 
is obtained by introducing a random phase variable. 
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I. INTRODUCTION 

Classical random walks (also known as 'drunken 
walks') have found practical applications in mathematics, 
physics and computational science, for example in stud- 
ies of diffusion, Wiener processes and search algorithms, 
respectively. Quantum physics introduces new perspec- 
tives, such as quantum diffusion quantum stochas- 
tics and quantum walks |^. The quantum 

walk is particularly appealing as an intuitively accessi- 
ble model underpinning quantum diffusion and quantum 
stochastics. Remarkable properties of these quantum 
walks (QWs) have been discovered; of particular inter- 
est is that the spread (standard deviation) for the quan- 
tum walk is proportional to elapsed time t, as opposed 
to \/t for the classical random walk; thus, the QW offers 
a quadratic gain over its classical counterpart. 

Physical implementations of the quantum walk have 
been proposed ||^, ^ , and possess the attractive property 
that they are inherently local in the sense that the spa- 
tial state shifts by one step along the lattice at each time 
step. One potential use of the QW is as a benchmark 
for assessing the non-classical performance of a quantum 
computer . ft is critically important in quantum infor- 
mation to develop algorithms and processes that behave 
in a distinct, observably different way than any classical 
one; the quantum walk is one such example. 

We extend studies of QWs to a higher number of spa- 
tial dimensions and examine the time dependence of the 
standard deviation, which reveals the universal feature 
of a quadratic gain over the classical random walk. We 
analyze and discuss the effects of entanglement between 
the different spatial degrees of freedom. We also compare 
with the equivalent classical random walk, and obtain the 
classical limit from the quantum model via the introduc- 
tion of a random phase variable at each time step and 
performing an ensemble average. 



II. THE ONE-DIMENSIONAL QW 

The classical random walk in one dimension describes 
a particle that moves in the positive or negative direc- 
tion according to the random outcome of some unbiased 
binary variable (e.g., a fair coin). The one-dimensional 



lattice on which the particle moves could be infinite or 
bounded (as in a circle). We may extend this to a QW 
by giving the particle an internal degree of freedom; for 
example, the particle may be a spin-1/2 system with 
internal Hilbert space 7^2 and basis states |±). The spa- 
tial state of the particle is given by a state in a Hilbert 
space Tispatiai of a one-dimensional regular lattice. Let 
with i an integer, denote the state of a particle lo- 
cated at position i\ the set {|i)} forms an orthonormal 
basis for Tispatiai ■ The total state of the particle is given 
by a state in the tensor product space 



'Ht — ^spatial ® ^2 ■ 



(1) 



Let the particle initially be in the spatial state |0) (i.e., 
localized at the origin) with internal state |-|-). To re- 
alize the 1-D QW 0], this particle is subjected to two 
alternating unitary transformations. The first step is the 
Hadamard transformation [p| , 



(2) 



which acts only on the internal state of the particle (i.e., 
on 7^2)^ and transforms the initial state |+) into the su- 
perposition --i= ( I -I-) -f I — ) ) . Following this transformation, 
we apply a unitary operator F that translates the posi- 
tion of the particle conditionally on the internal state: if 
the particle has internal state |-|-), it is moved one unit 
to the right, and if the internal state is |— ), it is moved 
to the left, i.e., 



r(N) 

F(N> 



(3) 



The translation does not alter the internal state, i.e., the 
states |±) are internal translation eigenstates. Since the 
transformation is linear, it will transform the superposi- 
tion state --i=(|-|-) -(- |— )) into a superposition state of the 
particle having moved left and right. Thus, the internal 
and spatial degrees of freedom become entangled. The 
Hadamard transformation is applied again, followed by 
F, and these transformations are repeated alternately. 

After n iterations, the particle is in an entangled state 
l^'n) G TLt- The probability Pi that the particle will be 
found at the i*'' location is given by 



P,. 



|((z|®(-|-|)|vI/„)|V|((z|®(-|)|vI;„) 



(4) 
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FIG. 1: The probability distribution of the 1-D quantum walk 
after 100 iterations. The internal state transformation used is 
the Hadamard transformation, and the initial internal state 
is h). 



In Fig. y, we plot the probability distribution of this 1- 
D QW as a function of i Analytical results are 

possible for the 1-D QW, and the n ^ oo asymptotic 
behaviour has been investigated Q . A key feature of the 
quantum walk is quantum interference, whereby two sep- 
arate paths between two nodes can interfere according to 
the phase difference. In contrast, the classical model has 
additive probabilities for alternate paths. Perhaps most 
interesting is the relative uniformity of the central por- 
tion of the distribution (—25 < i < 25) and the standard 
deviation of the distribution increases linearly with the 
number of steps t; this result is in contrast to the square 
root dependence of the classical random walk. Another 
peculiar feature is the asymmetry of the spatial probabil- 
ity distribution; this asymmetry is a consequence of the 
choice of initial state. The distribution resulting from the 
initial internal state \tfjs) = -75(1+) is symmetric. 



III. QW IN HIGHER DIMENSIONS 

The analysis of the one-dimensional walk can be ex- 
tended to higher dimensions. We define generalizations of 
the Hadamard gate, which place the internal state of the 
particle in superpositions of internal translation eigen- 
states, plus a generalization of F, which moves the parti- 
cle in the d-dimensional space conditional on the internal 
state of the particle. 

For a QW in d-dimensions, we require the particle to 
have an internal state in a 2''-dimensional Hilbert space. 
This internal state is simply described as the state of 
d coupled qubits thus, we can express the internal 



Hilbert space TYint 

Hint = ® «) • • • ® = ®'^n2 , (5) 

and give a basis for internal states in binary notation as 

\e1e2 ... Ed) = ® |e2) . . . ® led) , (6) 

where = ±. The state of the i'^ qubit (with basis 
|±)) will determine the direction (positive or negative) 
that the particle moves in the i*"^ dimension. That is, 
we define a translation operator F which translates the 
state of the particle by one unit in every dimension: the 
direction in the i^^ dimension is conditional on the state 
of the i*'' qubit. The internal translation eigenstates are 
those given in Eq. (^. 

For the 1-D QW, the quantum analogue of the classical 
"coin-flip" was the application of the Hadamard transfor- 
mation of Eq. (|^) . This transformation maps an internal 
translation eigenstate of the translation operator F (ei- 
ther |-f ) or |— )) into an equally weighted superposition of 
the two. The choice of phases in this transformation was 
to some extent arbitrary; the Hadamard transformation 
represents a choice with real entries. 

For the c?-dimensional QW, there exists a wide vari- 
ety of unitary transformations on the internal state that 
could be used as a generalization of the Hadamard trans- 
formation for the 1-D case. One obvious generalization 
would be to apply a Hadamard transformation H to each 
qubit in the decomposition of Eq. (^; i.e., the transfor- 
mation 

Hd = H (g) H (g) . . . (g) H . (7) 

This internal transformation is separable, in the sense 
that it does not produce entanglement between the spa- 
tial degrees of freedom. This choice could be viewed as 
the quantum analogue of using d independent coin tosses, 
one for each spatial dimension. 

Another obvious generalization, which is not separa- 
ble and does produce entanglement between spatial de- 
grees of freedom, is the 2'^-dimensional discrete Fourier 
transform (DFT) D^, defined as follows. Expressing 
the basis of Eq. (|^) as labelled by its numerical value 
/i = 0, 1, . . . , 2^* — 1}, the DFT acts on this basis as 

2''-l 

Note that the Hadamard transformation is the d — 1 dis- 
crete Fourier transform Di. As the Hadamard transfor- 
mation does for the 1-D case, this DFT transforms any 
internal translation eigenstate into an equally weighted 
superposition of all the eigenstates. Unlike the ten- 
sor product of Hadamard transformations, it is non- 
separable and highly entangles the different internal 
qubits. Although this internal transformation can also 
be viewed as a quantum analogue of d independent coin 
tosses, this entanglement between the spatial degrees of 
freedom is a genuinely quantum effect. 
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The DFT transformation represents a natural choice 
for the phase relationship between the translation eigen- 
states of the superpositions. However, this choice of 
phases is arbitrary, and we may consider other choices, 
which will have a different effect on the QW. We 
also investigate another internal state transformation 
(the Grover operator |^) that also produces an equally 
weighted superposition is the transformation (defined on 
the same basis as used above) 

G,|^) = ^(-2|^)+^|^)). (9) 

This choice, like the Hadamard transformation, possesses 
only real entries. 

There are, of course, an infinite variety of other non- 
separable choices for the internal transformation by em- 
ploying different phase relationships. Also, a bias could 
be introduced into the transformation, which would give 
an unequally weighted superposition of translation eigen- 
states; however, we consider only unbiased transforma- 
tions here. 

One of the remarkable properties of the 1-D QW is 
that, unlike its classical counterpart, it can produce an 
asymmetric distribution. Note, however, that with ap- 
propriate initial conditions (such as the state {ips) = 
-^(j-l-) -I- i|— ))) a symmetric distribution is obtained. It 
is of interest to question what effect the initial condi- 
tions will have on the higher-dimensional QWs. (Note 
that a symmetric distribution can always be obtained by 
averaging over initial conditions.) 

IV. CALCULATIONS OF QWS 

We begin our analysis with the straightforward gener- 
alization to higher dimensions of using the Hadamard 
transformation on each qubit. Fig. |l| and Fig. ^ 
show simulation results for the Hadamard walk both in 
one-dimension and a tensor product H (g) H for two- 
dimensions respectively. The initial condition for the 
internal state was chosen to be the separable state com- 
posed of all qubits in the |— ) state, which leads to an 
asymmetric probability distribution. 

For the case of separable transformations with separa- 
ble initial conditions, the different spatial dimensions be- 
have independently; thus, the variance can be expressed 
in terms of the one-dimensional case. For example, con- 
sider the family 

H, H®H, H(g)H(g)H, ... ; (10) 

the time dependence of the standard deviation for these 
walks is plotted in Fig. ^ and the corresponding slopes 
Acr/Ai are presented in Table |. We observe that 

Acti Aa2 A0-3 _ Aoi /^ Affi /^ Acri 
^ At' At ' At '■■■^"^ At At At 

(11) 



100 




^-foO -50 50 100 
X 

FIG. 2: Probability distribution of the 2-D quantum walk 
using the separable internal transformation H (g) H over 100 
iterations, with initial condition given by |— ) ® |— ). 
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FIG. 3: Time dependence of standard deviation for the series 
H, H (g) H and H ® H ® H with initial state given by ®''|-). 

where ad is the standard deviation for the d-dimensional 
QW, as expected for independent distributions. Also, 
by calculating a d-dimensional QW using this separa- 
ble internal transformation and projecting the state onto 
the Hilbert space for any one dimension, the state of 
the 1-D QW is recovered. Again, this property illus- 
trates that the different spatial dimensions are indepen- 
dent. Analytical results for these QWs follow from the 
1-D case in a straightforward manner. 

We now consider the behaviour of higher-dimensional 
QWs that possess entanglement between the spatial de- 
grees of freedom, i.e., QWs that have non-separable in- 
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Transformation Act/ At 



H 

H( 
H( 



H 

.H( 



H 



0.4544 ± 0.0012 0.4544 ± 0.0012 
0.6427 ± 0.0017 0.6427 ± 0.0017 
0.7871 ± 0.0021 0.7871 ± 0.0021 



40 



TABLE I: The slope of the standard deviation as a function 
of time for the family (H, H (g) H, H (g) H (g) H, . . . ). The 
slope Act/ At is found by linear regression of data points where 
t > 10 (such as to allow stabilization of irregularities caused 
by initial condition), cti refers to the 1-D case. 
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FIG. 4; Probability distribution for the quantum walk using 
the d = 2 DFT (D2) over 100 iterations, with initial condition 
given by |-) (g) |-). 



ternal transformations, such as the DFT of Eq. (^) . Fig. ^ 
shows the spatial probability distribution of the QW 
with internal transformation given by the d = 2 discrete 
Fourier transform D2; note that this distribution is dis- 
tinct from that of the H (g) H QW. In particular, it has 
the feature that the density of the distribution is signifi- 
cant near the origin, in constrast to the separable H 
QW which possesses only average density at the origin. 
Note also that it is asymmetric for the initial condition 
|— ) ® |— ); the asymmetry appears to be a general prop- 
erty of the higher-dimensional QWs as it is for the 1-D 
case. 

The time dependence of the standard deviations for 
the d = 1,2,3 DFT walks are plotted in Fig. | and 
summarized in Table ||. In contrast to the separable 
case, the trend observed in the DFT family is AcTd/At = 
^J{d-\- l)/2A(Ti/At. This trend is in agreement with the 
three calculations {d = 1,2,3). For the family of DFT 
walks, the standard deviation grows linearly with time, 
but the slope is less than that for the separable case (the 
tensor products of Hadamard transformations); see Ta- 
ble |l| This suggests that the entanglement between the 
spatial degrees of freedom serves to reduce the rate of 
spread. 
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FIG. 5: Time dependence of the standard deviation for the 
Dd DFT series with initial state given by ®'^\—). Details are 
given in Table 



Transformation Act/ At 



V(rf+ l)/2ACTi/At 



Di (H) 

D2 

D3 



0.4544 ± 0.0012 0.4544 ± 0.0012 
0.5569 ± 0.0006 0.5565 ± 0.0015 
0.6449 ± 0.0007 0.6426 ± 0.0017 



TABLE II: Slope of the standard deviation as a function of 
time, and comparison to the suggested pattern. Act/ At is the 
slope found by linear regression of data points where t > 10 
(such as to allow stabilization of irregularities caused by initial 
condition). 



Choosing different relative phases in the internal state 
transformation can lead to vastly different distribu- 
tions. Fig. ^ shows the results of using the internal 
transformation G of Eq. (^) . This distribution is chara- 
terized by its marked localization at the centre, as well as 
possessing peaks at the "maximum distance" attainable 
in the number of iterations (100 units from the origin). 

Note that the time dependence of variance for the non- 
separable 2-D transformations (D2, G2) are quite simi- 
lar, although the probability density functions are quite 
different in appearance. (See Table III.) The choice of 
initial condition does not appear to have a significant ef- 
fect on the time dependence of the standard deviation. 



V. OBTAINING THE CLASSICAL RANDOM 
WALK FROM THE QUANTUM MODEL 

A classical distribution can be obtained from the quan- 
tum model by introducing a random element into the 
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FIG. 6: Probability distribution of the 2-D quantum walk 
with internal transformation given by G (see Eq. ^) over 100 
iterations with initial condition |— ) ® |— ). 



Transformation Initial State Aa/At 



(H®H) 




I-) 


0.6427 ±0.0017 


Da 




1-) 


0.5569 ± 0.0006 


Da 


1+)® 


1+) 


0.5569 ± 0.0006 


Da 






0.6234 ± 0.0005 


Da 


IV'-) 




0.6009 ± 0.0006 


G (see Eq. |) 


1-)® 


I-) 


0.5418 ± 0.0020 


G 


1+)® 


1+) 


0.5418 ± 0.0020 


G 






0.5988 ± 0.0006 


G 


IV--) 




0.5440 ± 0.0008 



TABLE III: Slope of the standard deviation as a function of 
time for various 2-D transformations and initial conditions. 
Aa/At is the slope found by linear regression of data points 
where t > 10 (to allow stabilization of irregularities caused 
by initial condition). Here \4's) ~ -^{\+) ± iK)) is the state 
that produces the 1-D symmetric distribution, and \ip-) = 
-^(|±) ® 1—) — 1—) ® 1±)) is the entangled singlet state. 



R(/3) 








(12) 



where <Tz — ( o -i ) is the Pauli spin matrix, and /3 G 
[0, 2tt). We then consider a QW where the phase between 
the 1+) and |— ) states is randomly selected at each inter- 
val from a uniform prior distribution over [0, 2tt). Rather 
than applying the Hadamard transformation as the inter- 
nal transformation, we apply 



H(/3) =R(/3)HR(/3)"i = 



1 



-i/3 



(13) 



with a phase /3 chosen randomly from the set [0, 27r) at 
each time step. The resulting distribution has a stardard 
deviation comparable to that of the corresponding bino- 
mial distribution, but exhibits strong interference effects. 
By averaging oven many trials, the distribution rapidly 
converges to the binomial distribution. 

These results for the one-dimensional case can easily 
be generalized to higher dimensions. For the separable 
d-dimensional QW, the generalization is straightforward: 
one simply replaces each Hadamard transformation in the 
tensor product with a random H(/3) at each step. The 
separability ensures that the resulting walk is equivalent 
to the 1-D walk in each dimension. 

For non-separable internal transformations such as the 
DFT, a straightforward extension is to apply 



(R.i(/?i) 



Ri(/3i)-i) 
(14) 



where . . . , (3d are random phases, each from the set 
[0, 27r). That is, an independent random phase is added 
for each dimension (qubit). Again, by averaging 400 
walks of 50 iterations each, we obtained the 2-D bino- 
mial distribution to a high degree of confidence. 



VI. CONCLUSIONS AND DISCUSSION 



transformation at each time step. As shown previously, 
the "quantum" behaviour of the QW is due to the phase 
relationship (interference) between the separate paths of 
the walk. By adding a random element to the phase and 
averaging over many trials, we show that the quantum 
inteference can be made to disappear and that the dis- 
tribution of the classical random walk is regained. The 
introduction of this random phase is an example of de- 
coherence. 

Let us first investigate the one-dimensional case. In 
the internal translation eigenstate basis |±), the unitary 
operator that transforms the relative phase between these 



We present here a framework for calculating and an- 
alyzing quantum walks in higher dimensions. The gen- 
eralization of these walks beyond one dimension gives 
a wide variety of choice for the phases involved in the 
"quantum coin toss" . We discuss the role of entangle- 
ment between the diiferent spatial degrees of freedom as 
a possible non-classical property of the higher dimen- 
sional QWs. As different choices lead to different spatial 
probability distributions, it may be that specific unitary 
transformations of the internal Hilbert space are partic- 
ularly well suited for certain computational tasks. 

As with the one-dimensional QW, the increased rate 
of spread (given by the linear dependence of the standard 
deviation on time) is present in the higher dimensional 
walks. This property may be particularly valuable for 
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classical random walk based algorithms, such as quantum 
searches. We show that entanglement between the spatial 
degrees of freedom reduces the slope of this linear growth 
but not the linear dependence on t. These results are 
shown to be independent of the initial internal state in 
the cases investigated. 

We show that the classical distribution can be obtained 
from the QW by introducing an internal transformation 
with a random phase and then averaging over many tri- 
als. This result is expected; the quantum behaviour of 
the QW is due to interference effects between the phases 
of different paths. For higher dimensional QW, more 
random parameters (one for each spatial dimension) are 



needed. 
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